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ABSTRACT 

: 

We study nonlinear wave phenomena in self-gravitating fluid systems, with a partic- 
>- ■ ular emphasis on applications to molecular clouds. This paper presents analytical results 

for one spatial dimension. We show that a large class of physical systems can be described 
by theories with a "charge density" q(p); this quantity replaces the density on the right 
hand side of the Poisson equation for the gravitational potential. We use this formulation 
^ ' to prove general results about nonlinear wave motions in self-gravitating systems. We 

Q\ ■ show that in order for stationary waves to exist, the total charge (the integral of the 

charge density over the wave profile) must vanish. This "no-charge" property for solitary 
waves is related to the capability of a system to be stable to gravitational perturbations 
for arbitrarily long wavelengths. We find necessary and sufficient conditions on the charge 
density for the existence of solitary waves and stationary waves. We study nonlinear wave 
motions for Jeans- type theories [where q(p) = p — po] and find that nonlinear waves of 
>■ . large amplitude are confined to a rather narrow range of wavelengths. We also study 

wave motions for molecular clouds threaded by magnetic fields and show how the allowed 
range of wavelengths is affected by the field strength. Since the gravitational force in 



one spatial dimension does not fall off with distance, we consider two classes of models 
with more realistic gravity: Yukawa potentials and a pseudo-two-dimensional treatment. 
We study the allowed types of wave behavior for these models. Finally, we discuss the 
implications of this work for molecular cloud structure. We argue that molecular clouds 
can support a wide variety of wave motions and suggest that stationary waves (such as 
those considered in this paper) may have already been observed. 



Subject headings: hydromagnetics - wave motions - interstellar: molecules - stars: for- 
mation 
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1. INTRODUCTION 



Molecular clouds are self-gravitating fluid systems and are capable of supporting a 
wide variety of oscillatory motions. The overall goal of this work is to understand the 
dynamics of volume density waves in molecular clouds and to explore the extent to which 
waves can explain the observed cloud structure. These clouds comprise a substantial 
fraction of the gas in the galaxy, they provide the initial conditions for star formation, 
and they are interesting astrophysical objects in their own right. Molecular clouds are now 
fairly well-observed and exhibit very complicated highly nonlinear structures (see, e.g., 
the recent review of Blitz 1993). Thus far, however, no definitive theory exists to explain 
or predict the observed large scale structure of these clouds. Our working hypothesis 
is that volume density wave motions can produce (at least in part) this complicated 
structure. 

The study of nonlinear waves and solitons began over a century ago (see the classic 
papers by Russell 1844 and Riemann 1858). Similarly, the study of gravitational instabil- 
ity in astrophysical fluids has had a long and distinguished history (see Jeans 1928). In 
the meantime, however, remarkably little work has been done on the intersection of these 
two subjects, i.e., the study of nonlinear waves and solitons in systems where self-gravity 
is important. Liang (1979) has studied nonlinear waves in the cosmological fluid and 
has shown that soliton solutions are not allowed for ordinary pressure laws (see also Ray 
1983). Gotz (1988) has searched for solitons in Newtonian gravity, although the solutions 
found there are not relevant for our present discussion. Other related previous work does 
not directly address the problem of nonlinear wave motions in a self-gravitating fluid. 
Nonlinear Alfven waves have been studied and have been proposed as a mechanism for 
producing clumpy structure in clouds (Elmegreen 1990), although this study does not 
include self-gravity. Much of the previous work has been devoted to the linear stabil- 
ity of clouds (see, e.g., Langer 1978; Pudritz 1990; Dewar 1970; see also the review of 
Shu, Adams, & Lizano 1987). On a smaller size scale, the formation of molecular cloud 
cores has been considered through the process of ambipolar diffusion (Mouschovias 1976, 
1978; Shu 1983; Nakano 1985; Lizano & Shu 1989; Shu et al. 1987). On the size scale 
of molecular clouds, the study of wave motions including self-gravity remains largely 
untouched. 

Self-gravity is an important ingredient for the existence of stationary nonlinear waves 
in a neutral fluid. In the absence of self-gravity, acoustic waves in fluids are known to 
steepen and shock (see, e.g., Whitham 1974; Shu 1992). On the other hand, self-gravity 
provides dispersion (Jeans 1928) which tends to spread out wave packets. Thus, fluids 
with self-gravity can, in principle, reach a balance between nonlinear steepening and 
gravitational dispersion. This balance leads to the possibility of nonlinear stationary 
waves and solitary waves, which we study in this paper. 

In an earlier paper (Adams & Fatuzzo 1993; hereafter Paper I), we began a study of 
nonlinear waves and solitons in molecular clouds. We found one class of nonlinear waves 
for molecular clouds with no magnetic fields; we also showed that no soliton solutions 
exist for these systems. The one dimensional system of Paper I (with no magnetic fields) 
is highly idealized and hence unphysical in the following ways: (1) The usual Poisson 
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equation for gravity in one dimension produces a gravitational force that does not fall 
off at large distances, i.e., there is no asymptotic regime where the gravitational force 
goes to zero. (2) The one dimensional formulation of the problem does not allow for any 
two (or three) dimensional effects. (3) Magnetic fields are ignored, but are known to be 
dynamically important. (In Paper I, we briefly considered clouds with magnetic fields and 
derived a model equation which allows a wide variety of solutions, including nonlinear 
waves, solitons, and topological solitons - see §2.3 below). (4) Real (observed) molecular 
clouds are not collapsing as a whole, whereas the one dimensional model does collapse. 
In other words, the one dimensional model does not allow for a static equilibrium state. 
We note that items (1) and (2) are coupled in the sense that any physical system will 
have a finite spatial extent in the direction perpendicular to the wave motion. This finite 
extent implies a finite mass and hence allows gravity to fall off with distance. Similarly, 
items (3) and (4) are also coupled. Magnetic fields are (at least in part) responsible for 
the fact that clouds are not collapsing as a whole. 

This present paper is a generalization of Paper I. Since the full magnetohydrodynamic 
problem is rather complicated and the chances of finding meaningful analytic solutions 
seem remote, we use the full physical problem as motivation to derive a collection of 
model equations which approximate the true behavior of the system. In order to obtain 
analytic results, some of these approximations are by necessity rather severe. However, 
our goal is to obtain a qualitative physical understanding of the behavior of nonlinear 
waves in self-gravitating fluids. This approach has the advantage of allowing for analytic 
results which in turn provide us with a clear picture of "what depends on what" . Since 
we must sacrifice quantitative accuracy in our model equations, this work should be 
complemented by both numerical studies and by studies in higher spatial dimensions. 
We are currently pursuing these avenues of research. 

This paper is organized as follows. We begin with a general formulation of the 
problem in §2. We develop the concept of a "charge density" where we replace the 
right hand side of Poisson's equation for the gravitational field with a more complicated 
function (this concept was introduced in Adams, Fatuzzo, & Watkins 1993; hereafter 
AFW). We show that many types of physical systems can be modeled in this manner. 
In §3 we prove general results which hold for all theories with a charge density. In 
particular, we prove that stationary waves must have zero total charge when integrated 
over one wavelength. We also derive necessary and sufficient conditions on the charge 
density for the existence of stationary waves and solitary waves. In §4 we consider the 
special case where the charge density q = p — p ; this case corresponds to the traditional 
Jeans analysis. In §5, we include the effects of a magnetic field which is perpendicular 
to the wave propagation direction. In §6, we study nonlinear waves for theories where 
gravity is modeled with a Yukawa potential. This formulation of the problem allows the 
gravitational force to fall off at large distances. In §7, we show how stationary waves 
in two spatial dimensions can produce an effective charge density theory in one spatial 
dimension. We conclude in §8 with a discussion and summary of our results. 

2. GENERAL FORMULATION FOR ONE-DIMENSIONAL WAVES 

In this section, we introduce a class of model equations for the study of nonlinear 
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waves in molecular clouds. We begin with the equations of motion for a molecular cloud 
fluid. In one spatial dimension, these equations take the form: 



du du 1 dp dip 
dt dx p dx dx 

d 2 ib 

J = 4«G P . (2.3) 

For most of this paper, we take the pressure of the molecular cloud fluid to have a general 
barotropic form, 

p = p{p), (2.4a) 

which includes most equations of state of interest. In order to obtain explicit results, we 
sometimes must specify the form of the pressure law. In these cases, we adopt the form 

p = a 2 s p + po \og(p/p R ) , (2.4b) 

where the first term corresponds to an isothermal equation of state (a s is the isothermal 
sound speed). The second term arises from a "turbulent" contribution. This term is 
motivated by the observational finding that linewidths in molecular clouds vary with 
density according to Av ~ p~ x l 2 (see, e.g., Larson 1981; Myers 1983; Dame et al. 1986; 
Myers 1987; and Myers & Fuller 1992). If we interpret the observed linewidth as the 
effective transport speed in the fluid, we obtain (Av) 2 = v 2 uih = dp tUT b/dp = po/p, 
where the final equality follows from the observed relation. We then obtain the form 
given in equation [2.4b] by integration (see also Lizano & Shu 1989; Myers & Fuller 
1992). 

We find it convenient to work in dimensionless units. We let pr denote a refer- 
ence density and a s denote the sound speed. We can non-dimensionalize all quantities 
according to 

u — > u/a s , (2.5a) 
P -> p/pR, (2.5b) 
x — > kx where k 2 = AtiGpr/ a 2 , (2.5c) 
* -> ka s t, (2.5d) 

p - , (2.5e) 

^^^/a 2 s , (2.5f) 
• (2.5g) 

a 2 s PR 

For applications to molecular clouds, we are primarily interested in spatial size scales of 
1 - 30 pc where number densities are typically ~100 - 1000 cm -3 . The temperature 
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is expected to be in the range T = 10 - 35 K and hence the sound speed a s = 0.20 - 
0.35 km/s. The wavenumber k = 2it/Xj where Aj is the (usual) Jean's length that would 
result from thermal pressure (with sound speed a s and fiducial density pn) . For example, 
if we take a s = 0.2 km/s and pn = 2 x 10 -22 g cm -3 (the conditions roughly appropriate 
for the Taurus molecular cloud), we get x = 1 at a physical length scale of ~0.5 pc. The 
parameter k determines the relative size of the "turbulent" contribution to the pressure; 
the aforementioned observations indicate that k lies in the range 6 - 50. 

2.1 The Concept of Charge Density 

In this section we generalize the theory described in equations [2.1 - 2.4] by introduc- 
ing the concept of a "charge density" (see AFW). We can include a variety of additional 
behavior simply by modifying the Poisson equation to take the form 

where q(p) is a function of the density p. We denote the quantity q(p) as the charge 
density. We use this terminology because, in general, whatever appears on the right 
hand side of a Poisson equation is often called "the charge density". Notice that q(p) 
defined here has absolutely nothing to do with the electric charge density (which of 
course is completely negligible in this problem). Notice also that the choice q = p gives 
us back the usual system. In the following subsections (§2.2 - 2.4), we present physically 
motivated examples which can be written as charge density theories. 

One general motivation for introducing a charge density is to consider modifications 
of one-dimensional gravity (see below) and/or to include additional long range forces in 
the problem. We can model such systems by assuming the forces are conservative (~ 
VV'ex) and thus can be written in terms of a new potential ip ex . We then consider the 
potential ip to be the total potential, i.e., the sum of the usual gravitational potential 
and the new potential i/> ex . The dimensionless Poisson equation must have the form 

■Tr-t- = p + something = q , 
ox 2 

where we have defined the right hand side of the equation to be the charge density. Here 
we consider those special cases where the charge density is a function of p only. 



2.2 Jeans Theory 



The simplest possible nontrivial extension of the theory arises from subtracting out 
the contribution to the gravitational potential due to the background fluid (see, e.g., 
Jeans 1928; Binney & Tremaine 1987). This approximation can be written in the form 
of a charge density theory with 

q(p) = p-po, (2-7) 

where po is the dimensionless background density of the fluid. Thus, the traditional Jeans 
analysis results in a charge density theory. 
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One physically motivated way to obtain a theory with the same mathematical form 
as in equation [2.7] is to posit a physical system which rotates at a uniform rate O (see, 
e.g., Shu 1992). In the direction perpendicular to the rotation axis, Gauss's law for a 
uniform density cylinder with density po implies 

-rofi 2 = -2nGp w , (2.8) 

where w is the radial coordinate (note that all quantities in equation [2.8] have their usual 
dimensions). Thus, the uniform density state will be in mechanical balance provided that 

p = tt 2 /2nG. (2.9) 

It is straightforward to show (a homework problem in Shu 1992) that for waves propa- 
gating along the axis of rotation in this system, we obtain equations of motion with a 
"charge density" of the form q = p — po, where po is given by equation [2.9]. 

We now consider "typical" values for the rotationally induced fiducial density pq. At 
the solar circle, the rotation rate around the center of the galaxy is of order O ~ 10~ 15 
rad/s. With this value of O, the density implied by equation [2.9] corresponds to a number 
density of no ~ 1 cm -3 . In molecular clouds, the observed rotation rates are larger than 
this fiducial value by factors of 10 - 100 (see, e.g., Goldsmith & Arquilla 1985) and hence 
the implied number density becomes no ~ 10 2 — 10 4 cm -3 . These values thus lie in an 
interesting range for molecular clouds. 



2. 3 Charge Density from Magnetic Field Effects 

Another way to obtain a charge density is to simulate the effects of an embedded 
magnetic field. In Paper I, we derived a model equation on this basis. In this model, 
we assumed that the magnetic field points in a direction perpendicular to that of the 
wave motion and that the neutral component of the fluid is coupled to the magnetic field 
through its frictional interaction with the ionized fluid component (which is well-coupled 
to the field). We also ignored any chemical effects so that the ionized component of the 
fluid obeys a continuity equation. The resulting model equation can be derived from a 
theory with the charge density written in the form 

q(p) = p + T(--—), (2.10) 
\p p F J 

where V represents the coupling strength between the neutral and ionized components 
and where pp determines the ion density (assumed to be constant in Paper I). We note 
that the model equation corresponding to the charge density [2.10] is idealized in two 
ways. First, a dissipative term has been dropped to obtain this form. Second, the ion 
density is not calculated self-consistently. The net result of these approximations is to 
make the magnetic force on the neutral component into a long-range force that balances 
the long-range force of gravity. This treatment is unphysical in the sense that these 
magnetic forces, which arise from the frictional force exerted on the neutral components 
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by the ions, are intrinsically local. On the other hand, this approximation allows for 
something to cancel the long-range force of gravity 



2-4 Charge Density from Yukawa Potentials 

We can also derive a charge density for theories in which gravity is modeled by a 
Yukawa potential. The motivation for this approximation is to include the effects of a 
decreasing gravitational field strength while retaining the one-dimensional treatment of 
the problem. As we show here, Yukawa potentials provide one means of realizing this 
behavior. In particular, we write the Poisson equation in the generalized form 

-^=m 2 ^ + p. (2.11) 

The Green's function for the operator d 2 /dx 2 — m 2 has an exponential fall off and hence 
produces an exponential fall off in the gravitational force between point masses. The 
value of the parameter m determines the effective range of the force. ' 

The right hand side of equation [2.11] defines an effective charge density q, although 
it is not written explicitly as a function of density only. We note, however, that we can 
integrate the stationary version of the force equation [2.2] to obtain 

^ + h(p) + ^v 2 = E, (2.12) 

where E is the constant of integration (a discussion of the stationary wave approximation 
is given in §2.5). The quantity h(p) is the enthalpy and is defined by 



Hp) = fj- (2-13) 

The speed v can be eliminated by using the solution of the continuity equation for sta- 
tionary waves, i.e., v = A/p (see equation [2.17] below). If we now solve equation [2.12] 
for ij) and use the result in equation [2.11], we can read off the appropriate form of the 
charge density, i.e., 

q ( p ) = p + m *[E-h(p)-^]. (2.14) 

As we show in §7, another way to obtain a one dimensional theory where gravity falls 
off with distance is to begin with a two dimensional theory and then reduce it to a one 
dimensional theory through the use of rather severe approximations. We can thus obtain 
yet another charge density theory (compare equations [7.9] and [2.14]). 



' For completeness, we note that a more complicated additional term of the form m 2 ifj n 
could be used; this term represents nonlinear interactions in the gravitational field and 
is not of interest here. Notice, however, that this nonlinear theory can also be written as 
a charge density theory. 
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2.5 Stationary Waves and Solitary Waves 



In this paper, we are interested in stationary nonlinear waves and their limiting forms 
known as solitary waves and solitons. Since the definitions of these wave entities vary 
greatly in the literature (e.g., Whitham 1974; Coleman 1985; Rajaraman 1987; Drazin & 
Johnson 1989; and Infeld & Rowlands 1990), we present the following working definitions 
from Paper I: We use the term stationary wave to mean any wave structure that is a 
function of the variable £ = x — v$t only. We use the term solitary wave to refer to 
any solution of a nonlinear wave equation (or system of equations) which (1) represents 
a stationary wave, and (2) is localized in space so that the wave form either decays or 
approaches a constant at spatial infinity. We use the term soliton to be synonymous with 
solitary wave, although the term soliton is often reserved for waves which also satisfy 
the additional requirement: (3) the wave form can interact strongly with other solitons 
and retain its identity. Solutions which satisfy requirement (3) must have extraordinary 
stability in order to pass through each other and, after emerging from the collision, retain 
their initial forms. Such solitons are very rare; wave entities which satisfy the first two 
requirements and not (3) occur much more frequently. 

We first consider stationary waves, which correspond to traveling waves of permanent 
form. For these waves, the fluid fields are functions of the quantity 

£ >=x - VQ t, (2.15) 

where vq is the (nondimensional) speed of the wave. Next, we introduce a new velocity 
variable 

v = u — vq, (2-16) 

which is simply the speed of the fluid relative to the speed vq of the wave. Using the 
above definitions in the continuity equation [2.1], we find the relation 

pv = A = constant, (2-17) 

where the constant of integration A is the "Mach number" of the wave. 

We want to combine the equations of motion to obtain a single nonlinear differential 
equation for the density p. If we differentiate the force equation [2.2] with respect to x 
and use the generalized Poisson equation [2.6] to eliminate the potential, we obtain 



PPtt 



where subscripts denote differentiation and where we have eliminated the velocity depen- 
dence by using relation [2.17]. The equation of motion [2.18] is the fundamental equation 
of this paper. Notice its highly nonlinear nature. Fortunately and somewhat surprisingly, 
however, we can integrate this equation to obtain 



1 



pI = p 6 



-2 



P 2 ^-A 2 f{ P ) = Hn)- (2-19) 
dp 
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where we have defined f(p) to be an integral that depends on the form of the charge 
density q(p), i.e., 

The existence or non-existence of stationary waves for a particular physical system 
can be understood through the methods of phase plane analysis (see, e.g., Infeld & 
Rowlands 1990; Drazin & Johnson 1989). As illustrated above, in this method we reduce 
the system of equations to a single equation of the form 

^l = ^(p,C,), (2.21) 

where the Cj are constants (see equations [2.19] and [2.20]). Notice that equation [2.21] 
is just the virial theorem for an analogue particle moving in a potential — T (see, e.g., 
Rajaraman 1987); when viewed in this manner, equation [2.21] shows that the properties 
of the "potential" T determine the allowed behavior of the analogue particle and hence 
the properties of the wave solutions p(£) in a fairly simple manner. In particular, the 
form of J-{p) determines whether or not solitary wave solutions can exist. We also note 
that if dissipative terms are present in the original equation of motion, then the solution 
cannot be written in the form of equation [2.21] and hence stationary wave solutions do 
not exist. 

We first note that physically meaningful solutions must have T > (so that the 
solutions are real). For the systems considered here, the field p is a mass density and 
must always be positive. Thus, physically relevant solutions exist when T{p) is positive 
over a range of postive densities. Wave solutions exist when T{p) is positive between two 
zeroes of the function J 7 , where the zeroes of T correspond to maximum and minimum 
densities of the wave profile. The nature of the zeroes of T determines the nature of 
the wavelike solutions. For example, if T is positive between two simple zeroes, then 
(ordinary) nonlinear waves result. However, if T is positive between a simple zero and a 
double zero or higher order zero (i.e., any point where both T and dJ 7 /dp vanish), then 
a new type of solution - a solitary wave - can arise. Suppose we expand equation [2.21] 
about the double zero, which we denote as poo; we obtain 

p\ = {p- Po.?T" (Poo) + O [{p - Poo f] , (2.22) 

where T (p<x>) > because we are considering the case in which T is positive. Thus, as 
P —> Poo, the wave profile has the form 



p-Poo~*exp[=FV^'(Poo)£j, (2-23) 

where 5 is a constant. We see that a soliton consists of a single large hump of material 
and that the density smoothly approaches its asymptotic value p^ as £ — > ±oo. Formally, 
the wavelength of the solution diverges for a solitary wave. 

Another interesting type of behavior can arise when the function T is positive be- 
tween two double zeroes of J 7 , say pa and ps- In this case, the solution p(£) can approach 
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one value (e.g., pa) in the limit £ — > — oo and the other value (pb) in the limit £ — > oo. 
Solutions of this type are known as fcmfcs or topological solitons and are well studied in 
the context of quantum field theory (e.g., Coleman 1985; Rajaraman 1987). 

3. GENERAL RESULTS 

In this section we present general analytic results that apply to all theories of the form 
described in §2. These results are applicable for arbitrary forms of the charge density q(p). 
In particular, we state and prove four elementary "Theorems" which greatly constrain the 
allowed types of wave behavior for charge density theories. In the subsequent sections, 
we use these results to study the behavior of theories which contain specific forms for the 
charge density. 

3.1 The No-Charge Property for Stationary Waves 

To begin, we present an argument which shows that for any theory with a charge 
density q(p), a strong constraint must be met in order for physically relevant stationary 
waves to exist. This constraint arises from the fact that a stationary wave must have local 
extrema where the pressure and velocity gradients go to zero. Therefore, in order for the 
wave to remain stationary, the force of gravity must also vanish at these points. This 
behavior can only occur if the integral of the charge density q(p) between two extrema is 
zero. This argument can be stated as follows: 

Result 1. If a stationary wave solution p(£) exists for the one-dimensional theory, then 
the integral of the charge density over one wavelength must vanish, i.e., 

q = f x/2 de«[p(O] = 0, 

J-X/2 

where the wavelength of a soliton is taken to be infinite. This claim is limited to the case 
of nonsingular solutions; for solitons we also require that the density does not vanish at 
spatial infinity. 

In order to show that Result 1 is true, we first define the total charge Q contained in one 
wavelength to be 

,A/2 ,A/2 
Q= qdi = 2 qd£, (3.1) 

J-X/2 JO 

where A — > oo for soliton solutions. The second equality holds because of the symmetry of 
the problem about £ = (this result is valid for all classes of solutions except topological 
solitons, for which the formula for Q will be different by a factor of two). The first integral 
of the equation of motion for a stationary wave can be written in the form 

-l 

fl\ (3.2) 



Pi 



= ±V2p 



dp A' 
dp p 2 
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where the function / (defined in equation [2.20]) must be positive definite for valid solu- 
tions. Using the relation 



rf(/ 1/2 ) = \r 1,2q - 



a 2 dp 

p 2 dp 



dp, 



along with equation [3.2], we find the identity 

qd£ = ±V2d(f 1 / 2 ) 



(3.3) 



(3.4) 



Equation [3.4] may be substituted into equation [3.1], where care must be taken to ensure 
that the integrand has the proper sign. Using this result, we find 



rP2 

Q = 2V2 / d(f 1/2 ) , 
J pi 



(3.5) 



which can be integrated to obtain 



Q = 2V2 



f(P2) 1/2 ~ f( P l) 1/2 



= 0. 



(3.6) 



This result holds provided that we do not integrate over the singularity at the sonic point 
(see below). The second equality in equation [3.6] follows because the first integral of the 
equation of motion (and therefore /) must vanish at the values of density p\ (7^ 0) and 
P2 (by definition). Thus, the total charge Q must vanish. 

Result 1 greatly limits the allowed types of wave behavior in systems with self-gravity. 
For example, one important consequence of this "No-Charge Property" is that the charge 
density q(p) must vanish in the asymptotic limit £ — > 00 for a solitary wave or soliton. 
For these waves, the mass density approaches a constant as £ — > 00. 

In Paper I, we considered the case of molecular clouds with no magnetic fields and 
found solutions corresponding to stationary nonlinear waves. For this case, the charge 
density is given by q(p) = p and thus the total "charge" Q of the stationary wave is simply 
the total mass contained in one wavelength; this mass must be positive. Therefore, by 
Result 1, neither soliton nor stationary wave solutions can arise for this physical system. 

The apparent contradiction between Result 1 and the stationary wave solutions 
found in Paper I can be resolved by examination of their singularity. The stationary 
wave solutions of Paper I contain a singularity, p^ — > 00 as p — > pc, where pc satisfies 
dp/ dp = A 2 1 p 2 . Since v = A/p, the singularity is associated with the sonic point of the 
fluid, i.e., where the "effective" sound speed of the fluid is given by = dp/ dp. 

In the presence of the singularity, d(f 1 ^ 2 ) changes sign relative to dp across pc, and 
we must rewrite equation [3.5] as 



Q = 2^2 



PC 



d(f i/2 ) - 1 d(m 

Pi J PC 



P2 



•1/21 



(3.7) 
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The total charge is then given by 



Q = 2V2 



2/(pc) 1/2 -/(P2) 1/2 -/(pi) 1/2 



^0. (3.8) 



The discontinuity that allows these solutions to avoid the No-Charge Property necessarily 
produces a discontinuity in the gravitational force g(£), although the density profiles of 
the waves are continuous. As a result, the solutions of Paper I for the q(p) = p theory 
are unphysical. In this present work, we shall therefore require that realistic solutions of 
our equations be nonsingular. 

The most effective way to ensure nonsingular equations is to limit the discussion to 
a mass density range that excludes the sonic point pc- However, a nonsingular class of 
solutions can be found by choosing the constant of integration in equation [2.20] such 
that f(pc) = 0. For the theories considered in this paper, this class of solutions are 
characterized by a lower density bound of p\ = and a non-zero charge. We will consider 
this class of solutions briefly for the case of Jeans theory (see §4 and Appendix B). 

3.2 Relationship Between Solitary Waves and Jeans Stability 

In this section, we discuss the relationship between the existence of solitary waves and 
the absence of Jeans instability. This relationship arises because a soliton is in some sense 
an infinite wavelength perturbation on a uniform medium. On the other hand, the exis- 
tence of a Jeans length implies that all perturbations with a sufficiently large wavelength 
will collapse. Thus, the existence of a stationary perturbation of infinite wavelength (a 
soliton) is inconsistent with the presence of a Jeans length. This relationship can be 
stated more precisely as follows: 

Result 2. Suppose a physical system obeys a generalized model equation of motion of 
the form [2.18] and this model equation has solitary wave solutions. Then the system 
can have a uniform density state that is Jeans stable for arbitrarily large length scales. 

Suppose we have such a system. We can write the equations of motion in the form [2.1], 
[2.2], and [2.6]. We now consider a standard Jeans-type stability analysis, i.e., we take the 
unperturbed state to be p = po = constant and u = 0. (The existence of such a uniform 
density state requires q(po) = 0; systems which have solitary wave solutions always have 
such a zero - see the discussion of Result 3 below.) We expand the first order quantities 
according to 

f = fo + hj {kx -" t] , (3-9) 

where fi is a constant. After some algebra, the dispersion relation can be written in the 
form 

^(!)o^(|)„. < 3 - io > 

where the subscript denotes that the quantities are evaluated at the unperturbed density. 
Thus, the system can be stable to perturbations of all lengthscales provided that 

00 <0 (3.1D 
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for some density p = po > such that q(po) = 0. 

We now show that the condition [3.11] can be realized for systems which have solitary 
wave solutions. The generalized model equation of motion can be solved for q(p) to obtain 

q(p) = -e-^ 2 ^{e^(p)}, (3.12) 

where T is the first integral of the equation of motion (see equation [2.19]) and where p 
is an integrating factor, 

p = -61ogp + 2 log [p 2 ^ - A 2 ] . (3.13) 
We thus obtain an expression for dq/dp: 

|._ e -,/^ {e ^ )}+e -,^ {e ^ )}f |. (3 , 4 ) 

Since the system has solitary wave solutions (by hypothesis), we know that for those 
solutions the quantity e^jF has a double zero with a positive second derivative at some 
density p = po. Equation [3.14] shows that for p = po, the quantity dq/dp is negative. 
Thus, by equation [3.11], a Jeans stable configuration can arise for this physical system. 

We note that the converse of Result 2 is not true. Physical systems of this type can 
be Jeans stable to perturbations of all wavelengths and still not have soliton solutions. 
Thus, the requirement of Jeans stable configurations represents a necessary condition 
(and not a sufficient condition) for the existence of solitons. We consider other related 
conditions in the following section. 

3.3 Necessary and Sufficient Conditions 
for the Existence of Stationary Waves and Solitons 

In this section, we consider the conditions required for one dimensional systems to 
exhibit solitary waves and stationary waves. The question of whether or not solitary 
wave solutions exist is fundamental to the study of nonlinear dynamics and cannot, in 
general, be answered in a definitive manner. In this paper, we have shown that a fairly 
large class of physical systems can be modeled using a theory with a charge density q(p). 
The question can then be posed as follows: What properties must the charge density q(p) 
have in order for solitary wave solutions to exist and what properties are required for 
stationary waves to exist? Fortunately, we can provide a partial answer to this question. 

We first find a necessary condition on the charge density. In order for solitary wave 
solutions to exist, the first integral T of the equation of motion (see equation [2.21]) must 
be positive between a double zero and a ordinary zero (see §2). Using the solution in the 
form of equation [2.19], we can write 



■dp 



-2 



HP) = P P d P ~ A (3 - 15) 
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where f(p) is defined by equation [2.20] and depends on the form of the charge density. 
We restrict our discussion to waves which avoid the singularity at the sonic point. We 
also restrict our discussion to waves with finite velocity; this second restriction eliminates 
the point p = as a candidate for the double zero (if we let p — > 0, we get v = A/ p — > oo 
by the continuity equation). Thus, the required double zero of JF must occur at a double 
zero of f(p). Furthermore, f(p) must be a local minimum at the double zero. In addition, 
in order for the second (ordinary) zero of T to exist, the function f(p) must also have a 
local maximum (at some other density). By definition, the derivative of f(p) is given by 

d j_ = gM # _ * 

dp p p op 

Since we must avoid the singularity at the sonic point, both the required minimum and 
the maximum of f(p) must occur where the charge density q(p) vanishes. We denote 
the location of the minimum as p\ and the maximum as pu- In order for the point p\ 
to be a minimum, dq/dp must have the correct sign, which depends on the sign of the 
quantity in brackets in equation [3.16]. In other words, the sign depends on whether we 
are considering purely subsonic or purely supersonic waves. For subsonic waves we require 
dq/dp < at p\ and dq/dp > at pu- For supersonic waves, the sign requirements are 
reversed. These arguments can be summarized as the following necessary condition on 
the charge density q(p) for the existence of solitary waves: 

Result 3. In order for solitary wave solutions to exist for the class of theories considered 
in this paper, the charge density q(p) must have (at least) two zeroes p\ and pu- For 
subsonic waves, q(p) must be negative between the two zeroes; for supersonic waves, q(p) 
must be positive between the zeroes. (This result applies only to nonsingular solutions 
where pc does not lie between p\ and pu)- 

Result 3 is potentially very powerful as a test to see if solitary wave solutions exist. 
For example, many possible charge density functions do not have two zeroes and thus 
solitary wave behavior can be easily ruled out. We note, however, that Result 3 is not 
sufficient to guarantee the existence of solitary waves. We must place an additional 
constraint on the charge density to make sure that the second (ordinary) zero of the 
function f(p) exists. Given the two zeroes of q(p), the possible behavior of the function 
f(p) is sketched in Figure 1. The first possibility is that / has a second ordinary zero 
(solid curve in Figure 1) and thus solitary wave solutions exist. If the function / has 
another minimum at some density p% > pu and hence turns up (as shown by the short 
dashed curve in Figure 1), then we can always choose the constant of integration (3 to 
make p% a double zero of /. We thus obtain a depression soliton or void solution. Only 
for the special case where f(p) approaches a constant asymptotically (long dashed curve 
in Figure 1) does the system fail to have solitary wave solutions. 

For the usual case of purely subsonic waves, the existence of an ordinary second zero 
of / is guaranteed provided that the function / — > — oo as p ^ oo. In the limit p — > oo, 

df 1 dp 

Tp = - P dp q{ph (3 - 17) 
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where we have assumed that the equation of state is well behaved (i.e., the pressure 
increases with increasing density). Thus, the requirement that 



lim —£-q{p) > 

p^oo p Op 



(3.18) 



is sufficient, but not necessary, to ensure the existence of the second zero of f(p) for the 
case of subsonic waves. For most of the examples considered in this paper, q — > p and 
dp/ dp — > 1 in the limit p — > oo (for the case of a magnetic pressure with p ~ p 2 , dp/ dp 
— > oo); thus, the constraint [3.18] is always satisfied for these theories. 

We now consider the simpler case of ordinary stationary waves. The required condi- 
tions on the charge density q(p) for the existence of these waves can be stated as follows: 

Result 4. In order for stationary wave solutions to exist for the class of theories con- 
sidered in this paper, the charge density q(p) must have (at least) one zero pu- For 
subsonic waves, dq/dp > at pu &nd this zero must occur at a density larger than the 
sonic density pc (the density of singularity). Similarly, for supersonic waves, we must 
have dq/dp < at pu < Pc- 

The conditions outlined in Result 4 are both necessary and sufficient for the existence 
of stationary waves. This result is straightforward to understand. In order for such 
waves to exist, the function f(p) must have a maximum and two zeroes p\ < p2 such 
that pc < Pi for subsonic waves and pc > P2 for supersonic waves. The zero of q(p) is 
required to provide a critical point for f(p) - see equation [3.16]. The sign requirement 
on dq/dp ensures that the critical point is a maximum. The third condition, that the 
critical point pm be larger than the sonic point for subsonic waves and smaller than the 
sonic point for supersonic waves, is required so that the singularity does not lie in the 
range of densities in the wave profile. 



In this section, we explicitly consider the case of fluid systems with a charge density 
of the form q = p — p . As discussed above, this choice for the charge density is equivalent 
to adopting the original approximation of Jeans or to considering a uniformly rotating 
cloud (see §2.2). This system is also the simplest modification of the basic fluid equations 
that allows for the existence of stationary waves. Application of the No-Charge Property 
to this system requires po to be the average density for stationary wave solutions. 

We can analyze the possible behavior of this system using the results of §3. The 
charge density q(p) has only a single zero (at p = po) and hence Result 3 shows that no 
solitary wave solutions exist (except for solutions with vanishing density as £ — > 00; see 
Appendix B). In addition, dq/dp = 1 > everywhere and Result 4 implies that subsonic 
stationary nonlinear waves are possible provided that the singularity at the sonic point 
can be removed from the wave profile. This requirement implies that po > pc f° r subsonic 
waves and can be written as a constraint on the Mach number A: 



4. NONLINEAR WAVES AND SOLITONS IN JEANS THEORY 



A 2 



<Po 



2\ d P] 



(4.1) 
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This condition thus defines a maximum value of the Mach number A for nonlinear sta- 
tionary waves. 

In order to study the properties of the stationary waves, we must specify the equation 
of state and find the function f(p) defined by equation [2.20]. For the sake of definiteness, 
we take the pressure to have the form given by equation [2.4b]. For this equation of state, 
the integral in equation [2.20] can be evaluated to obtain 

f(p) =P- P -( K - P0 ) log(p) - + ^ , (4.2) 

where j3 is an integration constant. The function f(p) does allow for wave solutions; a 
schematic is shown in Figure 2. 

Without loss of generality, we can choose po = 1 (see Appendix A). This choice 
leaves us with a two parameter family of wave solutions. The parameter A sets the 
velocity of the wave and also determines the value of pc- The parameter A can vary 
from zero up to a maximum given by the condition [4.1]. For a given value of A, the 
constant of integration (3 merely slides f(p) in the vertical direction (see Figure 2) and 
thus determines the wave amplitude, which can vary from zero (when p\ = p2 = Po) to a 
maximum value (when p\ = pc)- In general, the range of possible amplitudes decreases 
with increasing A. A typical example of a nonlinear stationary wave in Jeans theory is 
shown in Figure 3. Note that the wavelength is A ~ 20 for this example; in physical units 
(scaled to the Taurus cloud - see the discussion following equation [2.5]), the wavelength 
is ~10 pc. 

Given the solutions for nonlinear stationary waves discussed above, we can determine 
the allowed wavelengths for nonlinear waves as a function of amplitude. For specified 
values of the constants A and /?, the wavelength of a stationary wave is given by 

A = 2 r (^y 1 dp, (4.3) 



pi 



an integral which unfortunately must be done numerically. In Figure 4 we plot the 
wavelength as a function of wave amplitude dp/p = (p% — Po)/ Po for two extreme cases. 
The first case has A = and corresponds to the "Jeans" length generalized for nonlinear 
waves; note that we recover the linear Jeans length for small amplitude. In fact, the Jeans 
length is shown to change very little over a large range of amplitudes. The opposite 
limit corresponds to the "shock" limit; the constraint that all velocities must remain 
below the effective sound speed a e ff gives a minimum possible wavelength for a wave of a 
given amplitude. This minimum wavelength approaches the "Jeans" length rapidly with 
increasing amplitude, so that for dp/p > 1 only a narrow range of wavelengths is possible. 

Nonlinear waves are characterized by relatively narrow, high density peaks. While 
the wavelength corresponds to the distance between peaks, the width of the peaks is also 
an interesting quantity. In order to obtain a measure of the width of the wave crests, we 
define the quantity X p to be the length of the region of the wave profile where p > po, 
i.e., 

" P2 Up 



Po 



x p = 2 h£ d P- ( 4 - 4 ) 
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In Figure 5, the width X p is plotted as a function of wave amplitude for both the "Jeans" 
and the "shock" limits described above. Figure 5 clearly shows that for 5p/p > 1, the 
range of allowed values of X p is very tightly constrained. In fact, the width X p has almost 
a single value as a function of dp/ p. 

These results have a relatively simple interpretation. For stationary waves, a balance 
must exist between dispersion, which acts to spread the wave, and nonlinear effects, which 
act to steepen the wave and cause it to shock. For the case at hand, the dispersion is 
supplied entirely by the self-gravitational interaction. (Recall that in the linear Jeans 
analysis the dispersion relation is given by uj 2 = a 2 s k 2 —AixGpo.) For waves with 5p/p > 1, 
nonlinear effects will be large. To balance these large nonlinear effects, the effects of 
self-gravity must also be large, hence the observation that the natural length scale for 
nonlinear stationary waves is of order the Jeans length. The narrow range of allowed 
stationary waves for large amplitudes 5p/p > 1 is thus a consequence of the fact that the 
amount of self-gravity needed to balance the nonlinearity is comparable to that required 
to cause the wave to collapse. 

For completeness we note that a class of solitary wave solutions exists for Jeans 
theory. In these solutions, the singularity is removed by choosing the integration constant 
(3 appropriately (see Appendix B for further discussion). An example of such a solitary 
wave is shown in Figure 6. 



5. NONLINEAR WAVES AND SOLITONS 
IN CLOUDS WITH MAGNETIC FIELDS 

In this section we consider the effects of magnetic fields, which are expected to play 
an important role in the dynamics of molecular clouds (see, e.g., Shu et al. 1987). In 
Paper I, we derived a model equation which incorporates the effects of magnetic fields in 
the equations of motion and produces a charge density theory (see the review in §2.3 and 
equation [2.10]). This model equation assumes that the largely neutral fluid is coupled 
to the magnetic fields through the ionic component as an intermediary. In order to 
obtain a simple form for the equation of motion, a dissipative term was dropped from 
consideration. Since this approximation has been discussed at some length in Paper I, it 
will not be considered here. 

Here we adopt an alternate approach; we assume that the magnetic field is strongly 
coupled to the neutral fluid. This assumption is valid for large lengthscales in molecular 
clouds (scales larger than those associated with dense cores), since the magnetic diffusion 
time is long compared to the magnetic sound crossing time (see Shu 1992). In fact, the 
existence of stationary wave solutions requires the absence of diffusive effects (recall that 
in Paper I we also ignored a diffusion term). 

With the neutrals strongly coupled to the magnetic field, the ions do not explicitly 
enter in the dynamics (the ions act only as an intermediate coupling agent). The equation 
of motion for the neutral component becomes 
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where we have considered only one spatial dimension and where B is written in dimen- 
sionless form B — > B /{4ita 2 pR) l t 2 ■ Observations of magnetic field strengths in molecular 
clouds (e.g., Goodman et al. 1989; Myers & Goodman 1988) indicate that B = 10 30 
fiG (with these values of field strength, the dimensionless version of B will be of order 
unity). The equation of motion [5.1] must be supplemented by an equation of motion for 
the magnetic field itself, 

— = V x (ux x B) . (5.2) 
dt v ' v ' 

Given the form of the magnetic force term, we see that the one-dimensional formalism 
is consistent as long as the magnetic tension does not result in a y-directed force (since 
no other forces exist to balance such a term). The magnetic field must also satisfy the 
usual divergenceless (no-monopole) condition 

V-B = 0. (5.3) 

These two requirements imply that the magnetic field cannot have components both 
parallel and perpendicular to x. Since a parallel component has no effect on the neutral 
dynamics, we consider a magnetic field of the form 

B = By. (5.4) 

Equations [5.1] and [5.2] may therefore be expressed as 

du du 1 d ( B 2 \ dip n 

m +u d-x + - p d- x V + ^) + ^ = ^ (5 - 5) 

and 

§ + A (uB) = a (5 . 6) 

In the stationary wave approximation, the equation of motion [5.6] for the magnetic 
field can be integrated to obtain the relation 

B v = constant . (5-7) 

This relation, along with the solution [2.17] to the continuity equation, gives the expected 
flux-freezing condition for a one-dimensional system, 

B . . 

— = a = constant . (5.8) 
P 

In typical molecular clouds, the magnetic contribution to the pressure is comparable to 
the "turbulent" contribution and hence we expect a 2 ~ k. We can use the flux freezing 
condition [5.8] in the equation of motion [5.5] to eliminate the magnetic field. The 
resulting equation of motion (including magnetic effects) is related to the previous one 
(equation [2.19]) through the transformation p — > p + a 2 p 2 /2. In other words, we have 
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simply added a magnetic pressure term to the nonmagnetic case. This new pressure is 
given by 

p(p) =p + Klog(p)+a 2 p 2 /2, (5.9) 

and we can proceed as in §4. For the values expected in molecular clouds, the second 
and third terms in equation [5.9] are comparable in size. In other words, the "turbulent" 
contribution to the pressure is comparable to the magnetic pressure. 

Here, we take the charge density q(p) to be the same as in the previous section. The 
results of §3 thus imply that the system cannot have solitary wave solutions, but can 
have subsonic stationary wave solutions. To study these waves, we find the function f(p) 
which becomes 

M =P-Y p2 + ^ p0 " !) ' " (" " po) Mp) - + |^ (5-10) 

where (3 is an integration constant. As we found in §4, the function f(p) is positive for 
small p and becomes negative as p — > oo. For this form of /, the wave solutions are 
analogous to those found in §4 and hence the discussion given there also applies to this 
present case. 

Figure 7 presents the density profile found by integrating equation [2.19] for the 
pressure given by equation [5.9]. The parameters used are representative of cloud envi- 
ronments. Comparison of Figure 7 with Figure 3 shows that a strong magnetic field is 
capable of sustaining a broader wave structure. This result is expected since the added 
magnetic pressure provides additional support against self-gravity and because the equa- 
tion of state [5.9] has a stiffer component. As in the case of pure Jeans theory (§4), 
nonlinear waves can only exist in a fairly narrow range of wavelengths. The allowed 
range of wavelengths for molecular clouds with magnetic pressure is shown in Figure 8. 
Notice that the maximum allowed wavelength varies more rapidly with amplitude than 
in the case of Figure 4. 

6. NONLINEAR WAVES AND SOLITONS WITH YUKAWA POTENTIALS 

In this section, we present a class of model equations with charge densities arising 
from modifications of the long range forces of gravity. In particular, we consider theories in 
which gravity is modeled with a Yukawa potential. As derived in §2.4, this approximation 
leads to a charge density theory with q(p) of the form 

q(p) = p - m 2 [h(p) + ^ - E] . (6.1) 
With this form for the charge density, the function f(p) can be easily found: 

M = P- p(p) ~Y + l m2 t h{p) + ^~ E ^- (6 - 2) 

Notice that for the equation of state [2.4b], the charge density q(p) is of cubic order and 
has either one or three real zeroes; thus, solitary wave behavior is possible. 
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We now show that all stationary waves in this theory must propagate at subsonic 
speeds. Result 4 shows that supersonic stationary waves must have q = and dq/dp < 
at some density pu < Pc- The first of these conditions (dq/dp < 0) implies that 



dp A 2 p 
dp p 2 m 2 



(6.3a) 



On the other hand, the requirement that the wave is supersonic (pm < Pc) implies that 



dp A 2 

IT' — <0 ' 
op p z 



(6.3b) 



where both conditions [6.3a] and [6.3b] are to be evaluated at the zero of q. Since both 
of these conditions cannot be simultaneously satisfied, no supersonic wave solutions exist 
for this theory. In other words, all physical stationary wave solutions must propagate at 
subsonic velocities. 

We now consider solitons in this theory. Result 3 shows that soliton solutions require 
the charge density q(p) to have at least two zeroes; for subsonic waves, q must be negative 
between two zeroes. Suppose q(p) has a zero at ps- We can always choose the constant 
of integration j3 to make the point ps a double zero of T . In this case, ps represents the 
fluid density far from the soliton. Without loss of generality, we can rescale the variables 
so that ps = 1 (see Appendix A). In order for the charge density to be negative between 
ps and a second zero of q, we must have 



dq 
dp 



1 - 



m 



Ps 



dp 
d~p 



A 2 ^ 
~72 



PS 



where the term in square brackets is to be evaluated at p 
condition for the existence of soliton solutions is that 



< 0, (6.4) 
ps = 1- Thus, a necessary 



m 2 > 



dp 
dp 



A' 



-l 



(6.5) 



This condition has two important consequences. First, we note that for the theory to 
have solitons at all, the quantity m 2 must satisfy the condition 



m 2 > 



' dp 




. dp 


i . 



(6.6) 



This constraint is equivalent to the condition that the theory has no Jeans length, as 
obtained from the linear dispersion relation (see also Result 2). Second, if the value of 
m 2 satisfies equation [6.6], then the Mach number A must satisfy the constraint 



A 2 < 



dp 


1 " 


[dp 


l m 2 



(6.7) 
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Since Aj p$ is the velocity of the soliton in the rest frame of the fluid, this constraint 
implies a maximum soliton velocity for a given value of to. 

The conditions outlined above are necessary for the existence of solitons. For equa- 
tions of state for which / — > — oo as p ^ oo, these conditions are also sufficient. Notice 
that this applies for our usual choice of equation of state [2.4b]. 

Thus far, our discussion has been limited to compression solitons, i.e., waves in which 
the density of the soliton is enhanced over that of the background. However, under certain 
conditions, depression solitons (voids) can exist as well. This type of wave behavior can 
occur if the function / drops below zero for pc < p < 1. For the equation of state [2.4b], 
we can show that a necessary and sufficient condition for the existence of depression 
solitons is that f(pc) < 0. After some straightforward algebra, this requirement can be 
written as a condition on the parameter to, 

2 ^ 2 [p(p c ) - p(l) - A 2 (l - p c ) 2 /2pl - h(p c ) + 

to < — o -. 6.8 

[h(pc)-h(l)+A*(l-f%)/2f%\ 2 

In Figure 9, we show wave profiles for both a standard (compression) soliton and a 
depression (or void) soliton. 

The constraints derived above define regions of allowed parameter space in the (to, A) 
plane. In Figure 10, we plot these constraints for the equation of state given in equation 
[2.4b] with k = 10. The vertical dashed line on the right of the figure corresponds to 
the sonic point; all wave solutions in this theory must be subsonic and hence A must lie 
to the left of this line. The constraints of equations [6.5] and [6.8] are shown as solid 
curves. In region I, ordinary (compression) soliton solutions are allowed; in region II, 
both compression and depression (void) solitons are allowed. If the constraint [6.5] is not 
satisfied, then no solitary waves exist (this case corresponds to region III of Figure 10). 
However, the theory does allow nonlinear waves of various amplitudes, which are deter- 
mined by varying f3. In analogy with the compression and depression solitons described 
above, these nonlinear waves can have profiles with narrow valleys as well as the usual 
profiles which are narrowly peaked. 

The inclusion of a magnetic pressure term a 2 p 2 /2 (see §5) can significantly change 
the results described above. For this case, h — > a 2 p in the limit p — > oo and thus 
/ — > (m 2 a 4 — a 2 )p 2 /2. For m 2 a 2 > 1, the function / — > oo as p — > oo in constrast to the 
behavior discussed above (see the discussion following equation [6.2]). For this pressure 
law, compression solitary wave solutions do not exist when m 2 a 2 > 1. 

7. APPROXIMATION OF TWO-DIMENSIONAL THEORY 

In this section we derive one dimensional model equations by approximating two di- 
mensional fluid systems. The goal of this procedure is to obtain a qualitative understand- 
ing of wave motions in higher dimensions while retaining the mathematical simplicity of 
the one dimensional theory. 

We begin with the equations of motion for a two dimensional fluid. We consider 
stationary wave solutions, so we define the variable £ = x — v^t and let v = u — vq be 
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the velocity along the x direction as before. We let w denote the velocity along the y 
direction. The continuity equation in this stationary wave approximation becomes 

ds(pv) + d y (pw) = 0. (7.1) 

The two components of the force equation become 

W£ + wv y + ti£ + ip£ = , (7.2) 

VW£ + WWy + hy + lj)y = , (7-3) 

and the Poisson equation for the gravitational potential is 

^« + $w = Qo(p) , (7-4) 

where we have left open the possibility of a nontrivial charge density qo(p)- 

We now consider a plane wave traveling in the x direction. However, we want the 
plane wave to have finite extent in the perpendicular (y) direction. As a first approxima- 
tion, we assume that the flow velocity in the transverse direction is much smaller than 
that in the direction of propagation and set w = 0. We are thus considering a wave 
structure with a more or less fixed profile in the transverse direction. In this case, the 
y component of the force equation [7.3] reduces to the hydrostatic condition h y = —i()y, 
and the Poisson equation can be written in the suggestive form 

V»« = Qo(p) + h w • ( 7 - 5 ) 



To make further progress, we must specify the equation of state (in order to specify 
the form of the enthalpy), the charge density qo, and the transverse profile of the wave. We 
note that the equation of state in the transverse direction need not be the same as that in 
the direction of wave propagation; in other words, the pressure need not be isotropic. For 
example, magnetic fields can provide support in the transverse direction. For simplicity, 
however, we consider an equation of state of the form [2.4b] and the enthalpy becomes 

h(p) = logp- n/p. (7.6) 

In order to isolate any possible pseudo- two-dimensional effects, we adopt a trivial form 
qo = p for the original charge density. If we assume that the wave structure is a localized 
lump in the transverse direction, we can take the density profile to have the separable 
form 

p = p(0 exp[-^ V] • (7-7) 

This form is chosen because it is symmetric about the y = plane and decreases with 
y. The parameter p, is like a wavenumber that sets the scale of the structure in the 
transverse direction. Using this ansatz for the density profile, we obtain a particularly 
simple form for the term h yy , namely, 

h yy = -2p 2 {l + K/p)-Ap A y 2 K/p. (7.8) 
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Because the equation of motion in nonlinear, we cannot self-consistently assume a profile 
of the form [7.7] for all values of y. However, near the center of the wave profile (i.e., near 
y = 0), we can neglect the y dependence everywhere except for keeping the first term in 
equation [7.8]. We thus obtain the Poisson equation in the form 

^ = p-2p 2 (l + K/p) = q(p), (7.9) 

where in the second equality we have defined the (total) effective charge density including 
two dimensional effects. We note that the form [7.9] is not sensitive to the particular 
form of the wave profile in the transverse direction; we simply require that the profile is 
even in y and decreases on a lengthscale given by p~ x . We stress that we have not solved 
the full two dimensional problem; we have simply assumed a reasonable structure in the 
transverse direction. The physical interpretation of this approximation is that if the wave 
structure has finite spatial extent in the transverse direction, then the gravitational field 
falls off with increasing distance (in the x direction). For the particular simple form of 
the transverse wave structure used here, we obtain yet another charge density theory, 
which can be solved using the methods developed in this paper. 

Given the form of the charge density (equation [7.9]), we see immediately that this 
theory does not allow for solitary wave solutions (because the charge density q(p) has 
only one zero with positive density - see Result 3). We also see immediately that dq/dp is 
always positive; thus, by Result 4, only subsonic nonlinear stationary waves are allowed. 
In order to ensure that physically relevant wave solutions exist, we still must show that 
the sonic point can be removed from the range of densities of the wave profile. This final 
requirement can be met if the zero of q(p) occurs at a density larger than that of the 
sonic point, i.e., provided that 

^ + p ^2 + 2k ] 1/2 > * + 1 ^ 2 + 4 ^ 2] 1/2 (7 _ 1Q) 

This condition can clearly be satisfied if the parameter p is large enough, i.e., if the 
density profile in the transverse direction falls off sufficiently rapidly. In order to obtain 
a better feeling for how restrictive this condition actually is, we consider the limit 4<k 
(typically, A ~ 1 and k ~ 10, so this limit is quite reasonable for molecular clouds). In 
this case, the constraint [7.10] reduces to the form 

where we have kept terms to leading order in A 2 /k 2 . Using the values A ~ 1 and k ~ 10, 
we find that p > 0.022 is sufficient to allow for stationary wave solutions to exist. In 
other words, the length scale of the density fall-off in the transverse direction can be 45 
times the thermal Jeans length and still allow for stationary waves (see equation [7.7]). 

Before leaving this section, we note that the Yukawa theory of the previous section 
can be recovered by using a pseudo-two-dimensional argument similiar to that given 
above. We write the Poisson equation as 

V>££ = P - Tpyy (7-12) 
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and then expand the equation about the plane y = 0. If we assume that the potential ip 
is separable near y = 0, we can write 

1>yy = VV>(£) • (7-13) 

Combining equations [7.12] and [7.13], we obtain the Poisson equation for a Yukawa 
theory (see equation [2.11]). 

8. DISCUSSION 

8. 1 Summary of Results 

In this paper we have developed further the theory of wave motions in self-gravitating 
astrophysical fluids. Although many of the results are general, the application to the fluid 
dynamics of molecular clouds is our primary motivation. Our results can be summarized 
as follows: 

[1] We have introduced the concept of "charge density" for the study of self-gravitating 
fluid systems (see also AFW). In this formalism, the density is replaced by a "charge 
density" q(p) on the right hand side of the Poisson equation (the continuity equation 
and the force equation remain as usual). We have shown that a large class of physical 
systems can be modeled with a charge density; we can thus include a wide range of 
physical effects while retaining a simple semi-analytic theory. 

[2] We have proven a "No- Charge Property" (Result 1) which shows that no solitons or 
stationary waves can exist in one-dimensional self-gravitating fluids unless the total 
charge vanishes (where the total charge is the integral of the charge density over one 
wavelength). This requirement greatly constrains the types of model equations that 
allow for stationary wave behavior. 

[3] We have shown that in order for a physical system to exhibit solitary wave behav- 
ior, the system must also be capable of having a configuration which is stable to 
gravitational perturbations of arbitrarily large wavelengths (see Result 2). We have 
thus discovered a fundamental relationship between gravitational stability and the 
existence of solitary waves. 

[4] We have found constraints on the form of the charge density q(p) for the existence 
of solitary waves and ordinary stationary waves (see §3.3, Results 3 and 4). These 
conditions allow us to determine (or at least constrain) the possible types of wave 
behavior for any charge density theory without having to solve the equations of mo- 
tion. 

[5] The original Jeans analysis provides the simplest nontrivial theory with a charge 
density and has q(p) = p — po- F° r this theory, we have studied the propagation of 
nonlinear waves in one dimension and found a class of wave solutions. In spite of this 
theoretical idealization, these waves are probably more physical than the nonlinear 
wave solutions found in Paper I. 

[6] Using the charge density formulation of the problem, we have performed a nonlinear 
Jeans analysis for molecular clouds. We find that the Jeans length is a slowly in- 
creasing function of the wave amplitude dp/ p. The wavelengths for stationary waves 
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must be less than this Jeans length. In addition, we find that waves with sufficiently 
small wavelengths tend to shock and dissipate. As the wave amplitude Sp/p becomes 
larger, nonlinear waves become confined to a rather narrow range of wavelengths. 
These clouds thus select out a particular lengthscale as a function of density contrast 
(see Figures 4 and 5). 

[7] The effects of magnetic fields can be incorporated into this theory in two concep- 
tually different ways. In Paper I, we derived a model equation which takes into 
account the relative drift between the neutral and ionic species; this model equation, 
which omits a dissipative term, allows for a wide variety of nonlinear wave behavior 
including solitary waves and topological solitons. In this paper, we have studied the 
other extreme - flux freezing. In this case, the magnetic field simply adds another 
component to the pressure. As in the case of Jeans theory (see item [6]), nonlinear 
waves are confined to a narrow range of allowed wavelengths (see Figure 8). 

[8] We have studied theories with Yukawa potentials which allow the graviational force 
to fall off with distance while retaining a one-dimensional formulation of the prob- 
lem. This theory contains both stationary waves and solitary waves. In addition, 
depression solitons or voids are allowed. 

[9] Using a two-dimensional theory as a starting point, we have derived a model equation 
which takes the form of a one-dimensional charge density theory. We can thus heuris- 
tically take into account two-dimensional effects while retaining a one-dimensional 
formulation. This particular model does not allow solitary wave solutions, but does 
allow stationary waves. 

We have thus introduced a new class of "charge density theories" for the study 
of wave motions in self-gravitating fluids. We have proven general results (Results 1 
- 4) which provide powerful constraints on the allowed types of wave behavior in these 
systems. Finally, we have argued that many physical systems can be modeled with charge 
density theories and we have studied several examples relevant to molecular clouds. 

8.2 Applications and Comparison with Observations 

Many of the results of this paper are general and can be applied to any self-gravitating 
fluid system (and other similar systems with long range forces; see AFW). However, 
one primary motivation for this work is to understand the formation of substructure in 
molecular clouds. These cloud systems are self-gravitating (see, e.g., the reviews of Blitz 
1993; Shu et al. 1987) and exhibit highly nonlinear structures (e.g., Myers 1991; Houlahan 
& Scalo 1992; Wood, Myers, & Daugherty 1993; Wiseman & Adams 1993; Blitz 1993; 
Adams 1992; de Geus, Bronfman, & Thaddeus 1990). In addition, these clouds exhibit 
a wide range of velocity structure as measured by line of sight velocity variations and by 
varying line-widths. These observations indicate cloud motions which are faster than the 
thermal sound speed but generally slower than the Alfven speed. In the models of this 
paper, the fluid speeds are determined by the Mach number A which is normalized to 
the thermal sound speed. Thus, for applications to molecular clouds, A should lie in the 
approximate range 1 < A < (1 + k) 1 / 2 ~ 3. 

The first obvious signature of wave motions in molecular clouds is periodic or nearly 
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periodic structures. Several examples can be found in the existing literature. For example, 
the star forming region NGC 6334 contains five regularly spaced clumps, as deduced 
from the 69 fim continuum map of the region (see McBreen et al. 1979); the inferred 
wavelength is ~3 pc. Such behavior is not rare. A survey of 23 globular filaments (defined 
to be filamentary dark clouds which exhibit condensations) shows that "the most striking 
similarity among all of the globular filaments is the regularity of their segmentation" 
(Schneider & Elmegreen 1979). This survey also shows that the characteristic spacing 
of the condensations along the filaments is about three times the width of the filament. 
Another example of possible wave motion is provided by the Taurus Molecular Cloud 
complex, which shows evidence for a periodicity in the velocity of 21cm observations of 
self-absorption (Shuter, Dickman, & Klatt 1987). These authors interpret their data as 
velocity waves with a peak to peak amplitude of ~3 km/s and a wavelength of 32 pc 
(see also Gomez de Castro & Pudritz 1992). To summarize, periodic structures with 
wavelengths in the range 3 - 30 pc are often found in molecular clouds. 

The solutions presented in this paper represent stationary waves in the molecular 
cloud fluid. Fortunately, a characteristic observable signature of these stationary waves 
is the relation given by the continuity equation: pv = A. If we observe a candidate wave 
train in a molecular cloud, then the observed line-center velocity (along the propaga- 
tion direction of the wave) should vary inversely with the density. Notice also that the 
linewidth is known to vary with density (see the discussion following equation [2.4]). Since 
the linewidths decrease with increasing density, the observed linewidths, when measured 
along the direction of a wavetrain, should also be anticorrelated with the density. Keep 
in mind, however, that this linewidth variation is a signature of the equation of state, 
whereas the aforementioned line-center velocity variation is a signature of stationary 
waves. 

One example of stationary wave behavior in an observed molecular cloud may be 
provided by the Lynds 204 complex (McCutcheon et al. 1986). This cloud is highly 
elongated and shows periodic structure in column density with a wavelength of ~ 1° 
(~3 pc for an assumed distance to the cloud of 160 pc). In addition, these authors find 
that the line of sight velocities (of CO line observations) along the filament are very well 
correlated with the density enhancements in the sense that the most massive parts of 
the filament have the smallest velocity displacement from the mean vq (see especially 
their Figures 4 and 5). Finally, the authors derive the relation Mj(vj — vo) = constant, 
where Mj and Vj are the mass and velocity of the jth segment of the filament. This 
relation looks suspiciously like the solution to the one dimensional continuity equation 
for a stationary wave (see equation [2.17]). 

We have also found that nonlinear stationary waves select out particular length- 
scales for their wavelengths (see, e.g., Figures 4, 5, and 8). The observed wavelength 
should thus be a calculable function of the various physical parameters involved (wave 
amplitude, sound speed, magnetic field strength, etc.). Although several parameters are 
necessary to completely specify the preferred wavelength (or range of wavelengths) , all of 
the parameters are physical quantities and can be determined, in principle, independently 
of the the observations of the wavelengths. 

Another observable quantity is the relative widths of the waves. The candidate 
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wavetrains in molecular cloud regions appear as chains of "blobs" which are (more or 
less) lined up and regularly spaced. The aspect ratio of these blobs (at a chosen density 
contour) is thus a well-defined observable quantity (this aspect ratio has a value of 3 ± 1 
for the sample of Schneider & Elmegreen 1979). We are currently studying stationary 
waves on filaments and can obtain a prediction for the relative widths of the wavetrains 
(Gehman et al. 1993) as a function of the physical parameters of the problem. 

In addition to forming substructure within molecular clouds, wave motions might 
be partly responsible for the formation of molecular clouds themselves. Consider, for 
example, a large scale nonlinear wavetrain (with a wavelength A of several pc) traveling 
through low density (n ~ 1) atomic gas. Atomic regions with sufficiently large column 
densities can become self-shielding and transform into regions containing molecular gas. 
After the wave has passed by, a wake of molecular material can remain. This process 
would thus leave behind a permanent record of the passage of the wave. We note that 
somewhat similar behavior forms clouds in the Earth's atmosphere, where waves set up 
regions of both low and high density; in the regions of high density, moisture condenses 
to form clouds and the so-called herring-bone cloud structure is produced (see, e.g., 
the introduction of Infeld & Rowlands 1990). We note that molecular cloud formation 
probably involves many different physical processes (see, e.g., the review of Elmegreen 
1987). However, the large scale waves discussed here may play an important role and 
should be studied further in this context. 

8. 3 Discussion and Directions for Future Work 

This paper (see also Paper I and AFW) provides a first step toward an understanding 
of wave motions and structure formation in molecular clouds. However, much more work 
remains to be done. Perhaps the key point that one should keep in mind is that complex 
fluid systems, such as molecular clouds, can wiggle around in many different ways. In 
this paper, we have concentrated on the study of stationary waves (travelling waves of 
permanent form) in one spatial dimension. Within this class of waves, a great variety 
of solutions exists. However, an even greater variety of wave motions can arise from the 
time dependent problem. These more general types of wave solutions should be studied. 

Another related topic is the actual generation of the wave motions. As discussed in 
Paper I, we expect that self-gravitating clouds will tend to collapse and excite a wide 
spectrum of wave motions (see, e.g., Arons & Max 1975). The self-gravity of the cloud can 
provide a more than adequate energy source for the waves. Any wave motions which get 
excited will generally either grow or disperse. However, the waves which live the longest 
will be the waves of permanent form - the stationary waves studied in this paper. One 
classic example of this type of scenario was studied by Zabusky & Kruskal (1965). They 
considered the time evolution of the Korteweg-deVries equation (a model equation which 
describes, among other things, surface waves on water; see Korteweg & de Vries 1895) 
and showed that initially simple cosine waves develop into a train of highly nonlinear 
pulses of permanent form; these pulses are the classic example of true (stable) solitons. 
Analogous calculations should be performed for the self-gravitating fluid systems of this 
paper and for more detailed models which more closely resemble real molecular clouds. 

Wave stability provides another avenue for future research. Even though the waves 
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considered here are stationary and thus possess a balance between gravitational disper- 
sion and nonlinear steepening, they may be subject to instabilities (see, e.g., Pego & 
Weinstein 1992). In the likely event that the waves are unstable, we must find the wave 
configurations with the longest lifetimes; these waves will be the ones which provide, in 
part, the observed structure of molecular clouds. 
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APPENDIX A: SCALING TRANSFORMATIONS 



In this Appendix, we consider a class of transformations which leave the fundamental 
equation of motion of this theory invariant. We begin with the equation itself: 



9 dp A 



dp 



dp 2 . 



+ A(p) = , 



(Al) 



which is written with an arbitrary charge density. We now consider transformations of 
the density p and the variable £ of the form 

p^Ap, (A2) 

£ - 7£ • (A3) 
The equation of motion [Al] remains invariant under this transformation provided that 

7 = A~ 1/2 , (A4) 

and the constituents of the problem scale according to 

A — > AA , (A5) 

p^Ap, (A6) 

q -> Aq . (A7) 

The Mach number A and the parameters in the pressure p can generally be rescaled as 
required. In order for the scaling of the charge density (equation [A7]) to be satisfied, 
the parameters of the charge density must also be rescaled appropriately, as we discuss 
below. 
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For Jeans theory, q = p — po, and the rescaling of equation [A7] can be met if 
Po — > Ap - In practice, we choose A to make po = 1 and then rescale the remaining 
parameters as described above. 

Both the Yukawa theory and the pseudo-two-dimensional theory of §7 can also be 
rescaled. We require 

m 2 -> Am 2 , (A8) 

for the former and 

p 2 -> Ap 2 , (A9) 

for the latter. These scalings are not unexpected since both m and p represent inverse 
lengthscales (very roughly they are the wavenumbers in the direction perpendicular to 
the propagation direction of the wave) and since lengthscales transform as in equations 
[A2] and [A3]. 



APPENDIX B: SOLITARY WAVE SOLUTIONS IN JEANS THEORY 

In this Appendix, we consider a class of solitary wave solutions in Jeans theory (see 
§4 in the text). For these solutions, the singularity is removed by choosing (3 such that 
f(pc) = 0. We first expand / (see equation [2.20]) about pc to obtain 

where the C subscript represents evaluation of the term in brackets at pc- We now 
expand equation [Bl] about the point pc and obtain 



^#(---c(4)ir 



+ 0[(p-p c )]. (B2) 

For the particular equation of state used here (see equation [2.4b]), the expansion reduces 
to a simpler form 

\p\ = P \ (po - Pc) Ppc + k]- 1 +0[(p- pc)] . (B3) 

Since the term in square brackets is nonzero, and since / > near pc, the quantity p| 
must be positive and nonsingular. 

For physically valid solutions, the density must lie in the range between pi = and 
P2 (where the peak density obeys the ordering p 2 > po > pc)- The first integral of the 
equation of motion (see equation [2.19]) may be expanded about p = p\ = to obtain 

= + (B4) 
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Thus, in the limit p — > p\ = 0, the wave profile has the form 

'A 2 



Po 



1/2 1 



lei • (B5) 



It is evident from equation [B5] that this class of solutions represents solitary waves (see 
the discussion in §2.6). The density profile for this case is shown in Figure 6. The 
dashed lines correspond to the sonic points. This solution corresponds to an unshocked 
flow with both supersonic and subsonic regimes and with densities both above and below 
the "background" density po. 

One complication of this solution is that the boundary conditions in the limits 
£ — > ±oo are unphysical because the velocity v ~ p~ l ~ £ diverges. Furthermore, 
the total charge Q also diverges, as shown by substituting the limiting form of f(p — > 
0) — > A 2 p /2p 2 into equation [3.8]. In order for this solution to be physical, appropriate 
boundary conditions must be applied at some finite distance. Suppose we interpret the 
initial unperturbed fluid to be at uniform density po and to have zero initial charge. 
The nonlinear structures that form out of this configuration must also have zero total 
charge. Thus, one "natural" choice of boundary conditions is to assume a symmetric 
density profile that also has zero total charge, i.e., we cut off the density profile at a finite 
lengthscale such that the total charge inside vanishes. As shown by equation [3.8], the 
Q = solitary wave is bounded by the sonic points (since f(pc) = 0). 

It is worth noting that charge densities q(p) can exist for which these solutions (i.e., 
solutions in which the singularity is removed in this manner) produce wave profiles which 
are bounded by a lower density p\ ^ 0, thereby allowing physical boundary conditions 
at spatial infinity. 
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FIGURE CAPTIONS 

Figure 1. Schematic of the function f(p) showing possible behavior of a charge density 
theory in which q(p) has at least two zeroes. The solid curve shows the usual case in 
which a second ordinary zero of / exists. The long dashed curve shows the case in which 
f(p) asymptotically approaches a constant and hence no soliton solutions are allowed. 
For the remaining case shown by the short dashed curve, depression solitons are allowed. 

Figure 2. Schematic of the function f(p) for Jeans theory. 

Figure 3. Wave profile p(£) for nonlinear waves in Jeans theory. This example has A = 1 
and k = 10. 

Figure 4. Allowed range of wavelengths for nonlinear waves in the Jeans theory. The 
upper curve shows the (nonlinear) Jeans wavelength A j plotted as a function of amplitude 
5p/p of the wave. The lower curve shows the "shock limit", i.e., waves with wavelengths 
smaller than this critical value will shock the dissipate. The wavelengths are expressed 
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in units of the linear Jeans length for the pressure given by equation [2.4b] with k = 10. 
Nonlinear waves in this theory must have wavelengths which lie between the two curves. 

Figure 5. Allowed range of wave profile widths X p for nonlinear waves in the Jeans 
theory. The upper curve shows the width for the longest possible wavelength (the Jeans 
wavelength) plotted as a function of amplitude Sp/p of the wave. The lower curve shows 
the width for the "shock limit" . All quantities are the same as in Figure 4. The width of 
the wave profiles for nonlinear waves in this theory are thus confined to a rather narrow 
range. 

Figure 6. Solitary wave solution for Jeans theory for A = 3 and k = 10. Vertical dashed 
lines show the location of the sonic points (see §4 and Appendix B). 

Figure 7. Wave profile for nonlinear waves in Jeans theory including the effects of 
magnetic fields as a pressure term. This example has A = 1, k = 10, and a = 3; the 
constant of integration (3 has been adjusted to make the amplitude comparable to the 
field-free case shown in Figure 3. Notice the enhanced width of the waves relative to 
those of Figure 3. 

Figure 8. Allowed range of wavelengths for nonlinear waves in the Jeans theory with 
magnetic fields included in the flux freezing approximation. The upper curve shows the 
(nonlinear) Jeans wavelength Aj plotted as a function of amplitude dp/p of the wave. 
The lower curve shows the "shock limit", i.e., waves with wavelengths smaller than this 
critical value will shock and dissipate. Here, k = 10, a = 3, and the wavelengths are 
normalized as in Figure 4. Nonlinear waves in this theory must have wavelengths which 
lie between the two curves. 

Figure 9. Wave profiles p(£) for a solitary wave (upper curve) and a "depression soliton" 
(lower curve) in Yukawa theory. Both examples have A = 1 and k = 10. Notice that the 
depression soliton has a much lower density constrast than the ordinary solitary wave. 

Figure 10. Allowed regions of the A — m plane for various types of stationary waves (for 
Yukawa theory and k = 10). Solitary wave solutions are allowed in both regions I and 
II. Depression solitons are limited to region II. Ordinary nonlinear stationary waves are 
allowed in region III. The vertical dashed curve on the right side of the figure corresponds 
to the sonic point for this equation of state; all allowed Mach numbers A must be less 
than this value. 
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